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Weak reaction freeze-out constraints on primordial magnetic fields 
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We explore constraints on the strength of the primordial magnetic field based upon the weak re- 
action freeze-out in the early universe. We find that limits on the strength of the magnetic field 
found in other works are recovered simply by examining the temperature at which the rate of weak 
reactions drops below the rate of universal expansion (T w < H). The temperature for which the n/p 
ratio at freeze-out leads to acceptable helium production implies limits on the magnetic field. This 
simplifies the application of magnetic fields to other cosmological variants of the standard big-bang. 
As an illustration we also consider effects of neutrino degeneracy on the allowed limits to the pri- 
mordial magnetic field. 

PACS numbers : 98.62.En, 95.30.Cq, 98.80. Cq 



I. INTRODUCTION 



It is by now widely recognized that cosmological mag- 
netic fields could be a significant factor on almost every 
scale relevant to the structure and evolution of the uni- 
verse. Recent measurements of intergalactic mag- 
netic fields have provided evidence that magnetic fields 
on the order of a few 10 _6 G are ubiquitous, e.g. in 
numerous galaxies, galactic halos, and clusters of galax- 
ies. Similar large magnetic field strengths have also been 
found in protogalactic clouds ||. However, the origin 
of these fields and the existence of galactic and inter- 
cluster magnetic fields are still an issue of debate ||-|. 
For a long time, it has been supposed that a turbulent 
dynamo mechanism [Q may be at work in which the 
magnetic fields might have arisen from the exponential 
amplification of small seed fields by hydrodynamic tur- 
bulence. However, recent detailed models || show that 
this dynamo model is inadequate. The fields that can 
be generated in this process are an order of magnitude 
weaker than what is actually observed. Goldschmidt and 
Rephaeli (To) also showed that Faraday rotation mea- 
surements could not be explained by a turbulent dynamo 
model. Furthermore, the existence || of damped Ly-a 
lines in QSO's indicates that primordial magnetic fields 
existed at the early times. 

Faraday rotation measure (RM) is defined by the ro- 
tation angle of polarized light in a magnetic field. It 
depends upon the strength and spatial extent of the mag- 
netic field, the density of the associated plasma, and 
on the wavelength of the observed radiation 0. There- 
fore, the observation of Faraday rotation in radio sources 
both inside and behind clusters gives information on the 
strength of the inter-galactic magnetic field [Q . For ex- 
ample, the detected magnetic field in clusters of galaxies 
has a typical magnitude of a few 10 -6 G and a coherence 
length scale of 10 ~ 100 kpc pi[. Another important 



estimate of the cosmological magnetic field has been es- 
tablished by Kosowsky and Loeb |l2|. They argue that 
the existence of primordial magnetic fields at the last 
scattering surface may induce a measurable Faraday ro- 
tation in the polarization of the cosmic microwave back- 
ground radiation. According to their results, it should 
be possible to detect the presence of a magnetic field at 
photon decoupling which would corresponds to a present 
day magnetic field as small as B = 10 _9 G. 

General effects of magnetic fields in astrophysical and 
cosmological processes have been investigated by several 
authors |T^-[T^|. In particular, the effect of a primor- 
dial magnetic field on big-bang nucleosynthesis was first 
studied by Greenstein |I3] and Matese and O'Connell 
|l5| . They argued that if primordial magnetic fields 
of sufficient strength existed in the early universe, then 
these could have had direct influences on both the expan- 
sion rate of the universe and the nuclear reaction rates 
[|l4] [18[ . These influences also could affect weak reaction 
freeze-out and hence the abundances of the light elements 
produced during the primordial big-bang nucleosynthe- 
sis epoch. 

The weak interaction rates themselves can also be af- 
fected by magnetic fields. This can influence significantly 
the rate of production of 4 He and other light elements. 
Recently, Cheng, et al. (l(| have considered the effects 
of a primordial magnetic field on big-bang nucleosynthe- 
sis. In that work they used weak interaction rates mod- 
ified by the presence of magnetic fields. They obtained 
constraints on the maximum strength of the primordial 
magnetic field in the framework of standard big-bang nu- 
cleosynthesis. According to their results, the allowed pri- 
mordial magnetic field intensity at the end of nucleosyn- 
thesis (T w 10 8 K) is about B < 2 x 10 9 -10 n G on scales 
greater than 10 4 cm but smaller than the event horizon 
at the big-bang nucleosynthesis epoch. Subsequently, 
Grasso and Rubinstein p7| obtained an upper limit of 
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B < 10 12 G if the coherence length Lq of the magnetic 
field at the end of primordial nucleosynthesis is in the 
range of 10 << Lq << 10 11 cm. In another work, Ker- 
nan, et al. obtained an upper limit of B < 1 X 10 11 G 
at T = 10 9 K. 

The main purpose of the present letter is to point 
out that the essential results of those papers can be un- 
derstood from the condition of weak reaction freeze-out 
alone and not on the later nucleosynthesis epoch. This is 
important for two reason: one is that the magnetic field 
strength is being constrained at an earlier time (1 sec 
rather than 100 sec); the other is that the constraint can 
be easily obtained without the necessity to run a full nu- 
cleosynthesis computation. This will allow easy parame- 
ter studies of effects of magnetic fields on a wide variety 
of cosmological models. As an illustration we consider a 
neutrino degenerate model here. 

During early epochs, particle species and fundamental 
interactions depart from equilibrium as the temperature 
of the universe decreases [^9[|o) . One of the typical de- 
partures is that of weak interaction decoupling which oc- 
curs when the temperature of the universe was about 1 
MeV. Weak-interaction freeze-out has important direct 
effects on the evolution of the universe. The synthesis 
of light elements depends sensitively upon the neutron 
to proton ratio n/p at the start of the nucleon reaction 
epoch (T < 10 9 K). However, the n/p ratio is deter- 
mined by the competition between the weak interaction 
rates and the expansion rate of the universe. A higher 
equilibrium n/p ratio at weak-interaction freeze-out will 
lead to a higher 4 He abundance during big-bang nucle- 
osynthesis. 

In this work we constrain the strength of magnetic 
fields at the weak reaction freeze-out epoch, without con- 
sidering how such fields could have formed. Our purpose 
is, therefore, to estimate the maximum allowed strength 
of the primordial magnetic field which is still consistent 
with weak-interaction freeze-out at around 1 MeV. In 
order to determine the maximum allowed strength of 
a magnetic field at this epoch, we here estimate the 
change in weak interaction rates in the presence of a 
magnetic field as well as the change in the Hubble ex- 
pansion rate from the presence of magnetic energy den- 
sity. In addition, we consider the possible effects of neu- 
trino degeneracy on the primordial magnetic field. Al- 
though we calculate the main effects on weak reactions, 
we have neglected the much smaller effects from some 
higher order interactions, radiative corrections, and other 
medium effects on the weak interactions. The electron 
mass will also be changed in strong magnetic fields, i. e. 
m e (B) = m e (B = 0)[1 + <D(a e ) + •••], with a- ~ 1/137 
being the electron fine structure constant |21j| , but for 
the magnetic fields of interest, other terms higher than 
0(a e ) are negligible. We employ a system of units in 
which h = ks = c = 1, except when specific units must 
be attached to a result. 



II. EXPANSION RATE AND WEAK REACTION 
RATES IN THE PRESENCE OF A MAGNETIC 
FIELD 

The properties of an electron in an external magnetic 
field have been extensively studied in a number of pa- 
pers [^l|-p4| . In brief, the energy states of the electron in 
a magnetic field are quantized and the properties of an 
electron arc modified accordingly. In order to investigate 
the properties of an electron in a magnetic field, we must 
first solve the Dirac equation in an external static and 
homogeneous magnetic field. We make the convenient 
choice of gauge for the vector potential in which a uni- 
form magnetic field B lies along the z-axis. We then ob- 
tain the electron wavefunctions in a magnetic field as was 
calculated by Johnson and Lippmann in detail |p3|| . The 
dispersion relation for an electron propagating through a 
magnetic field is |l6| 



e = \pI 



2eBn s ] 2 + m e K, 



(1) 



where n s — n + ^ — s z , (n s = 0, 1, ...), n is the princi- 
pal quantum number of the Landau level, s z — ±1/2 are 
the electron spins, e is the electron charge, p z is the elec- 
tron momentum along the z-axis, m e is the rest mass of 
the electron, and k is the anomalous magnetic moment 
for an electron in the ground state (n = 0, s z = 1/2). 
For relatively weak fields (i. e., B < 7.575 x 10 16 G), 
k = -^c-^z. while for stronger fields, k = f^Hn^) 2 
p2|. However, the analyses available in the literature |2l]] 
indicate that these higher order corrections have a very 
small effect at the magnetic field strengths and temper- 
atures of interest. Thus, we can ignore the anomalous 
magnetic moment term. Strictly speaking, all states of 
the neutron and proton are affected in the presence of 
a magnetic field. However, the effect is also smaller by 
a factor (m e /M p ) 2 . Hence, the modification of the pro- 
ton and neutron states by a magnetic field can also be 
neglected. 

The main modification of the electrons in a magnetic 
field comes from the available density of states for the 
electrons p4| . The electron state density in the absence 
of a magnetic field is 



d 3 p 
(2tt) 3 ' 

In a magnetic field this is replaced with 



n s =0 



eB 



dp z 



(2) 



(3) 



This modification will affect the thermodynamic vari- 
ables for the electrons as well as the fundamental inter- 
action rates. 

Now let us consider the effects of a magnetic field on 
the expansion rate. If the primordial magnetic field is 
spread over sufficiently small distances compared with 
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the event horizon, the geometry of the universe is not 
affected |l4|,^j|. Thus, a Robertson- Walker metric is still 
appropriate and the expansion rate of the universe can 
be simply described by the usual Friedmann equation, 



H=^ = J^Gp, 
Rdt V 3 p ' 



(4) 



where G is the gravitational constant and p is the total 
mass-energy density 

p = p n + p e + p v + p b + pB- (5) 

The subscripts j,e,v,b, and B denote the mass energy 
due to photons, electrons, neutrinos, baryons, and the 
magnetic field, respectively. Thus, the presence of a mag- 
netic field changes the expansion rate in two ways. One 
is by the associated magnetic energy density 



PB = 



8tt 



7 



(6) 



where 7 = B/B c and B c = m 2 /e = 4.414 x 10 13 G is 
the critical magnetic field at which quantized cyclotron 
states begin to exist. The other is the modification of 
the electron mass-energy density. Since the phase space 
of electrons in a magnetic field is modified, the electron 
energy density is given by 



Pe(B) = 2 



(27T) 



7 V [2 - <5„ s0 ] 
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1 
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Vl+27n 3 



1 + e eo+< ^ 



(1 + 2 7 n s ) 
m e /T e 



(7) 

where e = E/m e , a = m e /T e , <f> e = u e /T e , and T e is the 
temperature of the electrons. Since the electron chemical 
pot ential is small in the early universe, i.e., <fi e < 10 -9 
[BOj , we can ignore it in the calculations. 

The weak reaction rates in a magnetic field have been 
derived by several authors [p^|-|l8|. In thermal equilib- 
rium, the inter-conversion between neutrons and protons 
is possible through the weak reaction processes: 

n + v < — > p + e~ (8) 
n + e + < — ► p + v (9) 
n< >p + e~ + D (10) 

These weak processes set the n/p ratio in various as- 
trophysical processes such as big-bang nucleosynthesis 
p9|p0[ , neutron star cooling EfJ, etc. The reaction rate 
for each process can be calculated using the well-known 
V-A theory. The total weak reaction rates for the con- 
version of neutrons into protons in an external magnetic 
field is written as 



00 „ 
T n _ p (B) = 1 [2 - S n , ] / 



de- 



n.=0 



Wl+2~fn s 



v/e 2 - (1 + 2 7 n s ) 
{e-q) 2 e ta 



(e + Q) 



\ + e b(e-q)-i e 



(11) 



where 1/r = g\ (1 + 3a 2 )m^/47r 3 ~ 3.26 x 10^ 4 sec -1 , 
g v = 1.4146 x 10" 49 erg cm 3 , and a = g A /gv - -1.262. 
The parameters used in Eq. (11) are defined as q — 
Q/m e , b — m e /T u , and £ e = p, Vi! jT v . The quantity 
Q = M n — M p ~ 1.293 MeV is the neutron-proton 
mass difference. £ e is the electron neutrino degeneracy 
parameter which remains constant in the expanding uni- 
verse, fi„ e is the chemical potential of electron neutrino, 
and T„ is the neutrino temperature. Here, we neglect the 
polarization of the neutron source fig] . 

The inverse total reaction rate for the conversion of 
protons to neutrons can be obtained from detailed bal- 
ance. 



(12) 



In the limit of vanishing magnetic field (7 — * 0), T n ^ p (B) 
reduces to r„_ +p (i? = 0) The weak interaction rates, 
Eqs. (11) and (12), increase as the field strength 7 in- 
creases for fixed temperature (Ref. |l7j and see Fig. 1). 

From the above we see that magnetic fields affect nu- 
cleosynthesis in two competing ways, through the expan- 
sion rate of the universe H and the weak reaction rates 
r„,. The magnetic energy density ps in the total en- 
ergy density accelerates the expansion of the universe 
and therefore increases the primordial 4 He abundance. 
At the same time, however, the weak reaction rates in- 
crease in a magnetic field. This extends weak-reaction 
equilibrium to lower temperature which reduces the 4 He 
abundance. 



III. WEAK-REACTION FREEZE-OUT IN THE 
ABSENCE OF THE NEUTRINO DEGENERACY 

As long as the reaction rates exceed the expansion rate, 
chemical equilibrium is obtained. The equilibrium ratio 
n/p in the absence of the neutrino degeneracy is then 
given by pp0| 



P 



= exp(~Q/T). 



(13) 



This ratio is maintained only as long as the reaction rate 
per baryon T remains greater than the cosmic expan- 
sion rate H. Then, at some temperature To, the weak 
reactions decouple. The neutron to proton ratio n/p sub- 
sequently remains nearly frozen at its equilibrium value 
exp(—Q/T]j). The n/p ratio does, however, slowly de- 
crease due to occasional weak reactions. Eventually the 
ratio is dominantly affected by the beta decay of free neu- 
trons which continue to decrease this ratio until all of the 
neutrons are bound into nuclei. The freeze-out tempera- 
ture Td of the weak interactions is essentially determined 
by the condition T(Td) = H(Td) for which one obtains 
T D fn 0.75 MeV in the absence of magnetic fields 

In big-bang nucleosynthesis, the abundance oi 4 He can 
be easily estimated by assuming that almost all neutrons 
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are incorporated into A He. The resulting mass fraction 
of 4 -He will be approximately 



Y„ 



2(n/p) 



BBN 



1 + (n/p) 



BBN 



(14) 



where (n/p) bbn is approximately determined by the 
equilibrium value of the neutron to proton ratio n jp and 
the neutron decay factor prior to deuterium formation. 
However, in the presence of strong magnetic fields (i.e., 
7 ~ 1), the neutron decays more rapidly (about 15 % 
faster) than in the field-free case (27). The neutron de- 
cay factor will thus decrease and the allowed freeze-out 
temperature Td from the primordial 4 He constraint will 
increase. Using this we can approximately estimate the 
allowed values of (n/p) bbn in the presence of primordial 
magnetic fields during big-bang nucleosynthesis. For ex- 
ample, adopting Y p ~ 0.25, we find T D ~ 0.77 MeV. 

Now if the weak reaction freeze-out temperature Td is 
fixed by the value of primordial Y pi we can obtain a limit 
on the primordial magnetic field strength at that epoch. 
In order to constrain the strength of the magnetic field 
we again use the condition 



T(B D )=H(B D ) 1 



(15) 



for a given weak reaction freeze-out temperature To- In 
Eq. (15), Bd denotes the magnetic field strength when 
the weak reactions decouple. Therefore 7d = Bn/B c . 

Figure 1 shows the weak interaction rates from neu- 
trons converting to protons T n ^ p and the expansion rates 
H as a function of 7 = B / B c for the given temperatures 
T = 0.6,0.8,1.0 and 1.2 MeV respectively and £ e = 0. 
For a given temperature, the expansion rates are seri- 
ously affected by magnetic fields for 7 > 1 . On the other 
hand, large magnetic fields (7 > 10) significantly affect 
the weak reaction rates. For To > 0.75 MeV, we can 
find the intersection points of H and Y. These are the 
magnetic field strengths 7d at which the weak reactions 
freeze-out. As a result the freeze-out temperature in- 
creases as 7 increases. 



IV. WEAK-REACTION FREEZE-OUT IN THE 
PRESENCE OF THE NEUTRINO DEGENERACY 

As an illustration of how the above analysis is easily 
applied to other cosmological paradigms, we now con- 
sider the effect of a possible chemical potential for the 
electron neutrino, /i„ e . In this case the distribution func- 
tions of the neutrinos are different from those of the anti- 
neutrinos. Here, we choose to ignore the muon and tau 
(to m and m T >> 1 MeV) chemical potential because at 
the nucleosynthesis epoch they have already disappeared 
through annihilation and decay. 

By this time there is no experimental or direct theo- 
retical constraints for the magnitude or sign of neutrino 
asymmetries. However, there exist indirect constraints 



on the neutrino degeneracies obtained from primordial 
nucleosynthesis [p8| , p9| . Including neutrino degeneracy 
has two effects: (1) the energy density of neutrinos in- 
creases, the expansion rate of the universe thus increases; 
(2) weak reaction rates are modified by the change in the 
electron-neutrino distribution functions. Therefore the 
non-vanishing electron-neutrino degeneracy can directly 
effect the equilibrium n/p ratio at weak reaction freeze- 
out. 

With such a nonzero chemical potential of the electron 
neutrinos, the energy density of neutrinos is given by E8] 



Pu_ 

Pi 



21 ( 

T 



10 (ie 
7 \ 7T 



I (k 

7 V 7T 



(16) 



The equilibrium n/p ratio is also related to the £ e by 
n 



- = exp[-Q/T - £ e ]. 
P 



(17) 



Eventually, increasing £ e leads to a smaller value of n/p 
when the weak reaction rates freeze out and hence a 
smaller production of primordial helium. 

There have been various estimates for the constraint on 
the neutrino degeneracy from big-bang nucleosynthesis 
p8| , j30 31 1 . With the constraints imposed by large scale 
structure formation, Kang and Steigman obtained the 
limits on the electron neutrino degeneracy of —0.06 < 
£ e < 0.15 H). Recently Ref. @ has shown that big^ 
bang nucleosynthesis calculations agree with the primor- 
dial abundances of light elements inferred from the ob- 
servational data if the electron neutrino has a chemical 
potential due to lepton asymmetry. They obtain the pos- 
sible constraints of neutrino degeneracy, 0.003 < £ e < 
0.083 for 3.1 < 7710 < 5.5, where n w = 77 x 10 10 , n is 
the baryon-to-photon ratio n^/n^. Similarly, Kim et al. 
fl3l| obtained |£ e | — 1 with recent updated constraints on 
primordial light elements. 

In this work, therefore, we will adopt £ e = 0.15 to 
see the maximum effect of neutrino degeneracy on the 
constraint of primordial magnetic field. We calculate Eq. 
(15) in both cases of £ e = and £ e = 0.15. Figure 2 shows 
the freeze-out temperature Td versus 70 at the weak re- 
actions freeze-out. We can see that 7d goes to zero as Td 
approaches to 0.75 MeV and there is no jd for Td < 0.75 
MeV. In the case of £ e = 0„ we can estimate the freeze- 
out temperature Td < 0.76 MeV which satisfies the pri- 
mordial 4 He constraint, Y p < 0.245 |32]. Therefore, we 
obtain j D < 0.14 for T D < 0.76 MeV. In the case of 
£ e = 0.15, the neutrino degeneracy raises the freeze-out 
temperature by ATd ~ 0.025 MeV at 7 < 1. However, 
the n/p ratio for the positive £ e is suppressed, hence Y p is 
reduced. Therefore we can allow a higher weak freeze-out 
temperature Td < 0.82 MeV (corresponding to 7d < 0.8) 
and still satisfy the constraint, Y p < 0.245. For the more 
stringent constraint £ e < 0.083, we have Td < 0.78 MeV. 
Thus, we obtain 7d < 0.2 for T D < 0.78 MeV and it is 
possible to have j D <0.8 for T D <0.82 MeV. 
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V. CONCLUSIONS 

In this work, we have shown that weak-interaction 
freeze-out in the early universe is sufficient to provide 
a constraint on the strength of the primordial magnetic 
field. If magnetic fields existed in the early universe, 
in particular, during the weak reaction freeze-out, they 
could have influenced both the expansion rate of the uni- 
verse and the weak reaction rates. Therefore the weak 
reaction freeze-out temperature would be changed with 
respect to the field-free case. Subsequently, these influ- 
ences affect the abundances of light elements which are 
produced in this magnetized environment. A chemical 
potential in electron neutrinos can change the equilib- 
rium n/p ratio [Eq. (17)] as well as increase the expan- 
sion rate of the universe. We also consider the possible 
effects of neutrino degeneracy on the primordial magnetic 
field limits. 

In order to determine the maximum allowed strength 
of any primordial magnetic field, we use the simple con- 
dition, T(Bd) = H(Bd) for the given freeze-out tem- 
peratures Try- Since the observational requirement for 
the primordial 4 He abundance is that Y p < 0.245 [j32"t , 
the weak freeze-out temperature Tjj should be less than 
about 0.76 MeV for £ e = or 0.82 MeV for £ e = 0.15. 
We therefore find an upper limits from our results of 

B<6.2xl0 12 G (18) 

for the weak reaction freeze-out temperature Tr> = 0.76 
MeV (£ e = 0) and 

B<3.5xl0 13 G (19) 

for the weak reaction freeze-out temperature Try = 0.82 
MeV (£ e = 0.15). 

Since the universe has been a good conductor through 
most of its evolution, the cosmic magnetic field will con- 
serve magnetic flux as it evolves p5| . We therefore can 
obtain a simple relation 



B cx R- 2 oc T 2 . 



(20) 



If we assume that the cosmic magnetic field has continued 
to rescale according to Eq. (17), our results imply that 
the present cosmic magnetic field is less than 5.8 x 10 _7 G 
(£ e = 0) and 2.8 x 10~ 6 G (£ e = 0.15L The latter is larger 
than the limit determined in Ref. @ (3 x 10 _7 G) by a 
factor 10 but the former is nearly similar. 

In this work we have shown that a constraint on the 
strength of primordial magnetic field can be inferred by 
weak-interaction freeze-out in the early universe without 
full calculations of big-bang nucleosynthesis. Similarly, 
another important decoupling phenomenon in the early 
universe is the neutrino decoupling which occurred about 
T ~ 1.4 MeV. Neutrino decoupling is also very significant 
to the evolution of the early universe. Clearly, this neu- 
trino decoupling will provide a constraint on the strength 



of primordial magnetic field before big-bang nucleosyn- 
thesis. This will be the subject of a forthcoming paper 
@. 
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FIGURE CAPTIONS 



FIG. 1 The total weak interaction rate for converting neu- 
trons to protons T n ^ p and the expansion rates H 
as a function of 7 = B/B c are plotted for the values 
of freeze-out temperatures Try = 0.6,0.8, 1.0,, and 
1.2 MeV and £ e = as noted. 

FIG. 2 Freeze-out temperature Tr> versus the magnetic 
field strength jd at the weak interaction freeze- 
out. The solid line corresponds to £ e = and the 
dashed line represents for £ e = 0.15. 
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